Effect of environmental spins on Landau-Zener transitions 
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Landau-Zener (LZ) transitions of a two- level system (e. g., electronic spin in molecular magnets) 
coupled to one or many environmental spins (e. g., nuclear spins) are studied. For rather general 
interactions the LZ problem is reduced to that of a Landau-Zener grid. It is shown analytically that 
environmental spins initially in their ground state do not influence the staying probability P. This 
changes if they are prepared in a statistical ensemble. For a more specific model with environmental 
spins in a transverse field, LZ transitions are studied in the case of well-separated resonances in the 
LZ grid. The full evolution of the system is described as a succession of elementary transitions at 
avoided crossings and free evolution between them. If the environmental spins are strongly coupled 
to the central spin, their effect on P is weak. In other cases LZ transitions are strongly suppressed 
and P is decreasing very slowly with the sweep-rate parameter e oc 1/v, v being the energy sweep 
rate. 

PACS numbers: 03.65.-w, 75.50.Xx, 33.25.+k 



I. INTRODUCTION 



Spin tunneling in molecular magnets^^ (see Ref. for 
a review) is strongly affected by the dipole-dipole inter- 
action (DDI) between the spins S of magnetic molecules 
(henceforth electronic spins or simply spins) and by the 
interaction of the latter with nuclear spins I (see, e.g., 
Refs. iHQ)- Since the tunnel splitting A is typically 
rather small, DDI and interaction with nuclear spins cre- 
ate an energy bias W on electronic spins that can largely 
exceed A and thus bring them in or off tunneling reso- 
nance. As the bias due to the DDI depends on the instan- 
taneous configuration of all electronic spins, spin tunnel- 
ing and relaxation in molecular magnets is a complicated 
many-body phenomenon. Also coupling to nuclear spins 
renders the problem a many-body character, especially 
if the coupling to electronic spins is strong enough and 
thus tunneling transitions of the electronic spins strongly 
perturb the nuclear subsystem. 

There was a big excitement about the experimentally 
observed slow nonexponcntial relaxation at low temper- 
atures in molecular magnets^ and a great effort to un- 
derstand it theoretically using the concept of "tunneling 
window"i°Ji±2i±2 defined by \W\ < A. The role of nu- 
clear spins is believed to be in creating an effective time- 
dependent quasirandom bias on electronic spins that can 
bring them into the tunneling window and thus allow 
them to relax via tunneling. 1 - However, a similar effect 
arises via the DDI due to tunneling of other spins. 

Another kind of experiments on molecular magnets 
is the so-called Landau-Zener (LZ) experimen t 15 ' 16 in 
which the external magnetic field is time-linearly swept 
through the resonance and the magnetization change re- 
sulting from spin transitions is monitored. The asymp- 
totic staying probability in the standard LZ effect is given 



by 
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where v = W = const is the energy sweep rate (contribu- 
tions of the four authors are reviewed in Ref. \2m . Relat- 
ing P to the observed magnetization change allowed to 
extract the splitting A for the molecular magnet Fes^ 5 - 
The results were in accord with theoretical predictions 
that the dependence of A on a magnetic field along the 
hard axis shows periodic suppression of tunneling.— In 
addition, experimental results for the effective splitting 
A (e) have demonstrated the existence of an isotope 
effect originating from different nuclear spins of the in- 
vestigated Fes isotopes.— A remarkable property of the 
LZ effect with linear sweep is independence of the re- 
sult for P of the damping of quantum levels^ that can 
exceed A. Another remarkable feature is that in the fast- 
sweep limit e <C 1 the lowest-order result P = 1 — e is 
robust with respect to the effects of the interaction dis- 
cussed below. This allowed successful interpretation of 
the experiments at fast sweep in terms of the standard 
LZ effect, except for the isotope effect*^ 

In the LZ setting, all spins go through the tunneling 
window at some moment of time, in spite of the inter- 
nal bias due to the DDI and nuclear spins. However, 
the internal bias has a nontrivial time dependence that 
makes the total energy sweep time-nonlinear and drasti- 
cally changes the LZ effect. In the quantum- mechanical 
language, crossing of two energy levels in the one-particle 
LZ effect transforms into a series of crossings of many en- 
ergy levels in the many-body case (se e, e.g., Fig. 1 of Ref. 
124 Fig. 2 of Ref. [H or figures in Refs. [2ril[27l . where many- 
body LZ effect is described by random matrices) . The re- 
sult of multiple quantum transitions at all crossings is in 
most cases to increase the staying probability P for the 
initial quantum state, compared to the non-interacting 
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case. Especially strong increase of P (that is, the reduc- 
tion of the transition probability 1 — P) takes place for 
the ferromagnetic coupling between the electronic spins. 
DDI can be considered as a mixture of ferro- and antifer- 
romagnetic interactions, the former dominating for the 
crystals elongated in the direction along the easy axis for 
electronic spins, that is usually the case. 

Analytical and numerical solution of the many-body 
LZ problem is in general a formidable task. In particular, 
a full quantum-mechanical solution is prohibitive because 
of a huge Hilbert space, 2 N for TV spins 1/2. Quantum 
mean-field approximation (MFA) 24 ' 28 makes the prob- 
lem tractable numerically, since the size of the numerical 
problem is only N. A toy model of the many-body LZ 
effect, N spins all coupled with the same strength with 
each other,— also leads to a numerical problem of size 
N, that allows to accurately study the effect of many 
crossings and to test the validity of the MFA. Probabil- 
ity for all spins to remain in their initial states (that is 
not related to the magnetization change) in the LZ model 
with pairwise interaction between spins can be calculated 
analytically. 29 

Effect of realistic interactions such as DDI on the LZ 
effect could only be studied in the fast-sweep limit, where 
a perturbative analytical solution at order e 2 can be 
found^ This allowed to theoretically explain the ex- 
perimentally observed onset of the many-body regime at 
slower sweeps.— Up to now no way to theoretically study 
the effect of DDI on the LZ effect in molecular magnets 
at arbitrary sweep rate could be proposed. 

The interaction with nuclear spins strongly reduces the 
effect of DDI on the LZ effect that otherwise would be too 
strong. In Mni2, there are 12 Mn nuclear spins 7 = 5/2 
in each molecule, coupled to the electronic spin S via 
the contact hyperfine interaction. Since for nuclear spins 
this interaction is much stronger than other interactions 
such as the Zeeman interaction with the external field 
and the nuclear quadrupole interaction, one can neglect 
the latter. Then nuclear spins do not have any dynamics 
and simply act on the electronic spins as frozen-in ran- 
dom bias fields. The ensuing inhomogeneous broadening 
of tunneling resonances strongly reduces the influence of 
spin-spin interactions on the LZ effect i^ 5 . 

In fact, there are a lot of other atoms in magnetic 
molecules, mainly hydrogen atoms, that possess nuclear 
spins. These nuclear spins interact with the electronic 
spins via the nuclear dipole-dipole interaction (NDDI). 
As the nuclear magnetic moment is by a factor of 10 3 
smaller than the electronic magnetic moment, NDDI is 
much weaker than DDI. However, because of a large num- 
ber of nuclear spins (120 protons in a Fes molecule 3 ^) 
their cumulative influence on the electronic spins is ap- 
preciable. In contrast to Mn nuclear spins in Mni2, m ~ 
teraction of the remote nuclear spins with the electronic 
spins is not necessarily dominating, hence these nuclear 
spins should have their own dynamics and the field they 
produce on the electronic spins is not frozen in. 

For a rather general electronic-nuclear spin Hamilto- 



nian, a perturbative approach up to second order in the 
bare tunneling matrix element is used in Ref. EH] which 
does yield the dependence on the nuclear spins which are 
described as fluctuating fields. 14 Restricting to a simplicr 
model, a renormalized splitting obtained by an instanton 
technique is used to calculate the transition probability 
up to second order. The result^ coincides with that for 
the fluctuating fields, i.e., no isotope effect exists in this 
approach. A similar approach by eliminating the nuclear 
spins within a coherent spin path integral formalism for 
the electronic spins has been performed in Ref. [32| to 
predict suppression of macroscopic quantum coherence. 
Since this was done without a sweeping field, no predic- 
tions for the LZ probability can be deduced. A purely 
numerical solution of the problem has been presented in 
Ref. |M 

The aim of the present paper is to clarify the role of 
nuclear spins possessing their own dynamics on the basis 
of standard quantum mechanics with no usage of stochas- 
tic arguments or functional techniques. We will be using 
the same simplest model as Ref. l3li an electronic spin 
coupled to several nuclear spins, that also experience an 
effective field in the perpendicular direction and are ini- 
tially in the thermally equilibrium state. Our approach 
consists of two steps. First we prove analytically that the 
LZ problem for a rather general electronic-nuclear spin 
Hamiltonian can be exactly mapped to that of a Landau- 
Zener grid consisting of two families of parallel ascending 
and descending lines vs the energy bias W. This allows 
one to obtain some general predictions without specifying 
the explicit form of the nuclear spin interactions. Then, 
in the second step, we construct an analytical solution 
based on the standard LZ scattering matrix at every ele- 
mentary crossing of the system's energy levels including 
quantum-mechanical phases. 

Problems of this kind have been considered 
previously . 31 ' 34 ' 35 ' 36 Demkov and Oshcrov have shown 
that in the case of one level crossing several energy levels, 
the staying probability after all crossings is a product 
of all elementary staying probabilities.— Noteworthy is 
that this conclusion is valid even in the case of crossings 
that are not well separated with each other. In the case 
of several parallel lines crossing another set of parallel 
lines, Brundobler and Elser made a conjecture that the 
same result, multiplication of probabilities, holds for 
every ascending \meM- Refs. 38 : j39| present derivations of 
the Demkov-Osherov-Brundobler-Elser (DOBE) formula 
with different methods. On the other hand, the prob- 
ability of transition to upper ascending lines is exactly 
zero since it would require evolution in a negative time 
direction (the so-called no-go theore m 38 ' 40 ' 41 ). 

For the probabilities of transitions to lower ascending 
lines and to descending lines, no general analytical result 
exists. Since there are different "trajectories" through 
different intermediate levels on the way between the ini- 
tial and final levels, transition probabilities oscillate be- 
cause of quantum- mechanical interference (see, e.g., Ref. 
|42T ). As in most cases these oscillations should be aver- 
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aged out because of the distribution of system's parame- 
ters, an incoherent formalis m 31 i 35 neglecting nondiagonal 
density matrix coefficients has been used. 

What concerns the problem with nuclear spins that 
we are considering here, one of our results is that nuclear 
spins initially in their ground state does not influence the 
asymptotic staying probability P for the electronic spin. 
On the other hand, here is a substantial increase of P 
if the nuclear spins are initially in their excited states. 
In the latter case and more generally for nuclear spins 
at finite temperatures, the staying probability exhibits 
an isotope effect. However, the effect disappears again if 
nuclear spins are strongly coupled to the electronic spin, 
as well as for arbitrary coupling constant in the fast- 
sweep limit. 

The paper is organized as follows. Sec. HT1 introduces 
a general model for a central spin 1/2 coupled to envi- 
ronmental spins and undergoing a LZ transition. The 
Hamiltonian of environmental spins can be arbitrary, in- 
cluding their interaction with each other. It is shown 
how the Hamiltonian of the system can be diagonalized 
with respect to the environmental spins, that leads to 
the Landau-Zener grid. In Sec. IIIII a particular model is 
introduced, the model with one nuclear spin / in an arbi- 
trarily directed effective field. This model can be easily 
diagonalized by choosing an appropriate frame for the nu- 
clear spin having the z' axis along the total field. In Sec. 
HVI thc staying probability P on the LZ grid is calculated 
for well-separated resonances by multiplying LZ scatter- 
ing matrices of elementary avoided crossings. In Sec. |V| 
the model is generalized for N equivalent nuclear spins /. 
This model can be reduced to the model with one nuclear 
spin as the states of the nuclear subsystem can be classi- 
fied in the total nuclear spin I tot that is dynamically con- 
served. In this section the effect of quantum-mechanical 
phase oscillations is considered. Averaging over initial 
thermal distribution of nuclear spins and over different 
Jtot leads to only partial averaging out of the oscillations 
for the simulated system sizes. Here the more realis- 
tic incoherent approximation is introduced that leads to 
smooth dependences P(e). The concluding section con- 
tains a discussion of the two main cases: Mn 55 nuclear 
spins in the Mni2 molecular magnet and nuclear spins of 
the protons. 



II. GENERAL FORMULATION 

Tunneling of an electronic spin S of a magnetic 
molecule under the barrier can be considered in the ap- 
proximation of two resonant levels as dynamics of a pseu- 
dospin 1/2 coupled to N environmental spins, say, nu- 
clear spins /. The Hamiltonian has the form 

H{t) = ~W(t)a x - ^Aa x - a z V ({I}) + H nuc ({I}) , 

(2) 

where {1} = I 1; . . . , Ijv, and a z , a x are Pauli matrices. 
We assume a time-linear energy sweep W(t) = vt with 



v = const. The third term here describes the contact 
hyperfine coupling or dipole-dipole interaction between 
the electronic spin and nuclear spins. The fourth term is 
the nuclear Hamiltonian that can also contain interaction 
between nuclear spins. 

Let us introduce basis states 



l£;{roj}) = |0®|{mj}), e = ±i, my 



(i) 



-L 



rr> 



(i) 
I > ■ 



(3) 

, m,j , that form an eigenbasis of a z 



{mi} 
and 4°: 

vzmmj}) = £|£;{m 7 }), 

= rnf \€;{mj}). 

The Hamiltonian matrix 

H(t) = {(f;{mt}|fl-(t)|f;{mi})} 



(4) 



(5) 



can be represented in the following block form corre- 
sponding to the two electronic states: 



H(i) = 

where H+_ = H_ + = 
nuclear indices, and 



H__(t) H_+ 
H+_ H ++ (t) 



(0) 



(A/2) 1, I being a unit matrix in 



I ±± (f)= T -^(t)I T V + H nuc . 

Here the (21 + 1) N x (21 + 1) N matrices with respect to 
nuclear indices are defined by 

V={({m z }|t/|{m z })}, (7) 

etc. 

Since H±± is Hermitcan, it can be diagonalized by a 
unitary transformation matrix U± 



U±M±±(t)UZ 1 = B±(t), 



(8) 



U>±(£) being diagonal. Applying the block unitary trans- 
formation 



U 



u_ o 

O ILL 



(9) 



to Eq. ©, O being zero matrix, one obtains the trans- 
formed Hamiltonian matrix 



"w- (v „:;,). e°> 



where 



V' = -(A/2)U_U+\ V't = -(A/2)U+UI 1 . (11) 

In the new basis |£;fc), there are two sets of nuclear 
states for the electronic spin up and down, £ = ±1. Di- 
agonal elements of the diagonal matrices B± (t) have the 
form 



Dt,k(t) 



W{t) 



(12) 
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k = 1, . . . , (21 + 1) N , sorted so that d^ tk increases with 
k. Elements of V are 



Vkk' 



(A/2) (U_U^ 



'kk' 



(13) 



Avoided crossing of (27 + l) N ascending (£ = — 1) lines 
with (27 + 1) descending (£ = 1) lines forms the 
Landau-Zener grid shown in Fig. [1] in a particular case 
of one nuclear spin 7 = 1. In general, all crossings are 
avoided crossings with splittings Akk' = |2VJf fc ,| , as de- 
fined by Eq. (p~5|) . Crossing of |l;fe) with |— 1; fe') occurs 
&tD+, k (t) = D-^(t) or 



W(t) = W kk > = d+, fe - d_, * 



(14) 



If the system is prepared at t = — oo in the state 
|— 1, k) , then the probability P kk to remain in this state 
is given by the DOBE formula as the product of staying 
probabilities at all crossings ! 34 ' 37 ! 38 ' 39 



Pfefe = exp 



exp 



- e £(A fcfc ,/A) 

k' 

-El( u - U ") fefe 



(15) 



With the use of 

Ei(u-u-) fcfe ,r 

k' 

= E( U - U ?V ( u - u 



one obtains 



(U-U+ 1 U + iJZ 1 ) kk = {l} kk = l 



Pkk = e 



(16) 



(17) 



the same as Eq. (TTJ). This result shows that P kk does 
not depend on nuclear spins, i. e., there is no isotope 
effect. However, since we are interested in transitions of 
the electronic spin alone, the relevant quantity is 



p = y^^n k p kk , 



(18) 



kk' 



where n k are populations of the nuclear levels in the ini- 
tial state. This formula assumes that initially nuclear 
spins are in a state described by a diagonal density ma- 
trix, e.g., a thermal equilibrium state. With account of 
the no-go theore m 38 ' 40 ' 41 and Eq. (fTT)) , the result reduces 
to 



P 



E 

k>k' 



n k Pi 



kk' 



(19) 



If all nuclear spins are initially in the ground state 
(n k = Sk,i), there are no terms in the sum and the stan- 
dard LZ result is reproduced. Otherwise P exhibits an 



isotope effect and increases above e~ £ , i.e., coupling to 
nuclear spins is hampering spin transitions. For not too 
low temperatures, nuclear spins are equidistributed in the 
initial state, i. e., n k — (21 + 1)~ N . 

The effect of hampering spin transitions is especially 
strong at slow sweep. Indeed, at slow sweep in the case 
of a standard LZ effect, the system practically follows 
the lower adiabatic level, so that the probability P of a 
transition to the upper adiabatic level (i. e., of staying 
on the accending diabatic level) is exponentially small. 
For the model with nuclear spins, the starting level of 
the system is in general not the lowest level, because 
nuclear spins are thermally distributed over their energy 
levels. In addition, splittings A kk / are widely distributed, 
so that there are very small splittings and the adiabatic 
limit is practically never reached. As a result, there are a 
lot of transitions between ascending and descending levels 
in both directions, so that the electronic spin performs a 
complicated motion on the Landau-Zener grid. 

Consideration in this section suggests that interaction 
between N nuclear spins does not change the situation 
qualitatively. Indeed, Eq. (TIT)]) has the same form with 
and without interaction, so that the topology of the LZ 
grid is not affected. Numerical results of Ref. [H show 
only a moderate effect of interaction between nuclear 
spins. Anyway, all feasible types of interactions between 
nuclear spins, such as the direct dipole-dipole interac- 
tion and indirect interactions via the electronic spins, are 
much weaker than the interaction of nuclear spins with 
the electronic spin and interaction of nuclear spins with 
the external magnetic field or with the gradients of the 
microscopic electric field^ 7 - due to the quadrupole mo- 
ment of nuclei. 

In the next section we will consider our basic model 
of an electronic spin interacting with one nuclear spin in 
an effective field that can be easily diagonalized. Gener- 
alization for the case of many nuclear spins will be done 
later in Sec. IVl 



III. THE MODEL HAMILTONIAN AND ITS 
DIAGONALIZATION 

Consider the Hamiltonian for an electronic spin cou- 
pled to a single environmental spin 



H 



1 , 



:W(t)a z 



1 



Aa x - Aa z I z - k z I z - A X I X (20) 



that is a particular form of Eq. ((2|). Terms with A'a z I x , 
A"a z I y , and A y I y can be added to the Hamiltonian but 
such terms can be eliminated by choosing its own sys- 
tem of axes x',y',z' for the nuclear spin. A stands for 
the contact hyperfine interaction or for NDDI between 
the electronic spin and protons, as said above. In Eq. 
([20)1 A z and A^ are energy-dimensional components of the 
field acting on nuclear spins, e.g., an external magnetic 
field, A a — g n [J, n H a , where a = x,y, fi n is the nuclear 
magneton, and g n is the nuclear Lande factor. A Mn 55 
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nucleus in the molecular magnet Mni2 has 7 = 5/2 and 
the magnetic moment A< Mn 55 = 3.45/i„, so that g n Mn 55 ~ 
3.45/7 ~ 1.38. A proton has 7 = 1/2 and the magnetic 
moment fi — 2.79^„, so that g H:P ~ 2.79/7 ~ 5.58. For 
the contact hyperfine interaction with the Mn 55 nuclear 
spins in Mn 12 one has A/k^ — 0.02 K. This interaction 
is very strong and equivalent to a magnetic field of 40 T 
applied to nuclear spins. 

For the description of the states of the nuclear spin it is 
convenient to use the basis of its eigenstates correspond- 
ing to the states |J.) and |f) , or £ = ±, of the electronic 
spin. Fixing the electronic spin in the state £ creates an 
effective field £A z e z on the nuclear spin then gives the 
effective nuclear Hamiltonian 



77^,Ti,cff — — • I, 



where 



F 5 = (£A + A z ) e z 



^X^X 



(21) 



(22) 



is the total field acting on the nuclear spin. This Hamil- 
tonian can be diagonalized by choosing the z'^ axis for 
the nuclear spin in the direction of Fg, i.e., 



£A + A Z A x 
H e, 

F 5 Fj 



(23) 



where 



F ( = ^{ZA + A z f + kl. (24) 
The vectors e z ^ are rotated away from e z by the angles 

£A + A z 



1% 



Ft 



(25) 



For A x <c A, A z and A z < A, (3, is close to and (3_ is 
close to 7T. For A x <C A, A z and A < A z , both /3 + and /3_ 
are close to 0. For the transverse nuclear axes we define 



£A + A Z 



F 



Ax 
Ft 



(26) 



The corresponding basis set of states has the form 

\Z,m I ) = \Z)\m I ) e , (27) 

where the rotated states \mi)f depend on £. Elements of 
the Hamiltonian matrix H' in the new basis, Eq. (fTOf . 
are defined by 



H' 



(28) 



With the help of I —I x > e x > + I y e y + I z > e z > one obtains 



Hi 



^ : m,, m ;^' (29) 



where 



E, 



-fW(t) - F smi 



(30) 



are the diagonal elements of H>± of Eq. 



A, 



£ {mi |m/)_ e 



(31) 



are the elements of V in Eq. (fTT|) for £ = — 1, and the 
elements of V'^ for £ = 1 . Elements of the Hamiltonian 
matrix are placed in order of £ changing from —1 to 1 
and mi changing from 7 to —I. 

The projector £ (mj can be expressed through 

the spin rotation matrix^ 



d ( XW = (m'\e~^\m) 



(I — m)\(I — m')\ 



(I + m)!(J + m')\ 



1/2 



(3 



m-\-m 



E 



{-l) n {I + m + n)\ 
(I — m— n)\(m + n — m')!n! 



sin — 

2 



(32) 



summation going over max(0, m! — m) < n < I — m. The 
final expression is good for f3 > 0. In the case /3 < one 

should use the relation d$ m (J3) = (-l) m '- m d# m (-/?). 
For large 7 numerical implementation of Eq. (|3^|) leads 
to precision problems. In this case it is much more con- 
venient to obtain d^) m (/3) numerically by finding eigen- 
states of the operator I z cos/3 + I x sin/3 and projecting 
them on \m) , eigenstates of I z . 



Using \mi}^ 



-Vih \mi) , where |mj) are the states 
of the initial basis, quantized along the z-axis, one ob- 
tains 



= -0t + /3„ e 

For A z 



| m'j) 



from Eq. (|25|) one obtains 
.4 



/3 = 2£ arcsin 



F' 



(33) 



(34) 



where F = A 2 + A x , and, in Eq. (|3l 



A 



(35) 



In Eq. (f2"9")l F 4iTOj describes 27+1 ascending (£ = -1) 
and 27+1 descending (£ = 1) lines. All crossings are 

avoided crossings, because all m , ^ 0. The sum rule 



E [€L(p) 



= i 



m'=—I 



leads to Eq. (fTT|) in our particular case. 

For A z = 0, one obtains the matrix 
(£ = —1, £ = 1) in the form 



V' 



A ( Af 4 
~2 



F F 



(36) 

in Eq. ((TOD 

(37) 
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FIG. 1: Landau- Zener grid for one nuclear spin 1=1 and 
A z = 0: (a) A > A*; (b) A < A*; 




7 / 



FIG. 2: Squares of splittings, (A' mj ) 2 a |d^> mj (/3)J on 

the LZ grid for a large "nuclear spin" / = 30 and f3 — n/4 
(A X /A = 2A1). 



relation between and A in Eq. (|20[). For A x ^%~> A 
the matrix V becomes nearly diagonal, that results from 
/?£ w /3_ ? thus /3 w 0. In the opposite case A x < A ma- 
trix V' becomes nearly antidiagonal, that results from 



VI 



7r. For A x yf> A level splittings 
are large between the levels hav- 
ing the same mj, other splittings are small, as shown in 
Fig. [TJ (b) . This means that transitions of the electronic 
spin mostly leave nuclear spins at their initial levels, a 
natural result for a small coupling to spins. For A x <C A, 
the maximal splittings are those between the levels with 
nuclear indices related by m'j = —mi, while other split- 
tings are small, as shown in Fig. [TJ (a). In the original 
basis with the quantization axis z, this means that the 
states of nuclear spins do not change, too, similarly to 
the other limiting case. Indeed, for A^ = there are no 
transitions of nuclear spins since the latter have no dy- 
namics. For A z 7^ 0, matrix V is cumbersome but has 
the same properties: Nearly diagonal for j3 w (small A) 
and nearly antidiagonal for f3 s» tt (large A) . 

For large I the distribution of splittings in the plane 
(to/, m'j) , defined by the spin rotation matrix d^ t (/?) , 
becomes more interesting, see Fig. O The splittings be- 
come very small outside an ellipse and maximal on the 
ellipse. Inside the ellipse, , (/?) ~ 1 but changes 

abruptly if the indices change by one. For = ir/2 the 
ellipse becomes a circle of radius /. For j3 — * or /3 — ► tt 
the ellipse degenerates into a vertical (m'j — mj) or hor- 
izontal (m'j = —mi) lines, respectively. For large /, even 
if (3 is close to (A x /A > 1) or it (A x /A < 1), there 
are many crossings inside the ellipse where splittings are 
comparable to A. 

In particular, for ft < 1 the matrix d^) m (P) of Eq. ([32]) 
becomes nearly diagonal and can be expanded directly as 



P (m' \il y \ to) 
P 



— Sjii'in ^ ( $ m ' , m + 1 ^ m , m + 1 
— 8rn' ,m — l^m,m — l) j 



(39) 



where l m ,m+i = yHJ + 1) m(m + 1). This is the case 
realized for A^ ^> A, where from Eq. (j3~i|) with £ = — 1 
follows /3 = -2A/F ^ -2A/A X and in Eq. (f^T]) 

^ — Vmim'j T" ^ ^mj.mj-l t mj.m/-l 

) . (40) 



for / = 1/2, 




in accordance with Eqs. ([3T)l and (f3"B")) in this limit. Sim- 
ilarly, the case A x <C ^4 corresponds to /? close to tt, 
so that d^, m ((3) becomes nearly anti-diagonal. The el- 
(38) ements adjacent to the antidiagonal are proportional to 
(A x /A) lmi,mi+x- I n both cases, diagonal or anti-diagonal 
form of d!£l (ft) requires IA/A X < 1 or IA X /A < 1, re- 



fer 1=1, etc. The matrix elements depend only on the spectively. 
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IV. TRANSITIONS ON THE LANDAU-ZENER 
GRID 

As in the general case considered in Sec. HJ the energy 
levels described by Eq. ([50)1 separate into two families, 
ascending lines for £ = — 1 and descending lines for £ = 
+1. There are crossings between these families at 



E + , +) (t)=E_ >m() (t). 



(41) 



Each ascending line crosses with 21+1 descending lines 
and vice versa. In the case A z = that will be considered 
below, crossings are defined by 



W 



k = - m { j +) = -21, -21+1,..., 21. (42) 
There are total 41+1 crossings. For a given k, the values 



of m\ ' and rrij satisfy 



m T — rrii ~ 



--/ •_ ra\ \m ( I +) < I. (43) 



dt hi j — k , 

For the leftmost crossing with k — —21 there is 

only one solution, {^m' I ' =— I,m^ = Ij , i.e., only 

one pair of lines is crossing here. Similarly, for the 
rightmost crossing with k — 21 there is only one 

pair of crossing lines with (ra\ "* = I,rrij = —Ij . For 



the crossing with k — —21+1, there are two pairs 
of crossing lines with |m| ^ = —I + l,m^ = I^j and 

(m\ ' = —I, TOj + ' = I — lj . For k = there are I cross- 
ing lines with all possible values of mj and the corre- 
sponding m^ +S> — mi- In general, the allowed values of 
m\ ^ and rrij are in the intervals 

max (-1 - k, -I) < m ( j +) < min(7 - k, I). (44) 



and 



max (-1 + k, -I) <m\^< min(7 + k, I) (45) 



In the case A <C A, A x crossings are well separated 
from each other. In this case the process consists of tun- 
neling transitions at the crossings and free evolution with 
phase accumulation in the ranges between them. Tran- 
sition at each crossing is described by the LZ scattering 
matrix 



M(A) 



P 



sign(A) v / T T7 7V 
-sign(A) v / l _r Pe^ VP 



(46) 

where P given by Eq. (JTJ) is the Landau-Zener staying 
probability and 

<P = tt/4 + ArgT (1 - iS) + 5 (InS - 1) (47) 

with 5 = e/ (2tt) is the scattering phase and T is Gamma 
function. 



Evolution of the wave function between level cross- 
ings reduces to the accumulation of the phase factors 
exp [«3>£, mj (t)] , where the phases are given by 



dt'E e , mi (t r ). 



(48) 



The change of the phase between the fcth and (k + l)th 
crossings is given by 



0+1/2) 



or, finally, 



(fe+l/2) 



1 

h 

1 F 

h v 



eF 2 

7rA^ 



ptk+1 


^vt + Fin i 


/ dt' 







F (2k + 1) + Fmi 



2mj 



(49) 



(50) 



Evolution of the state on the interval between the kth 
and (k + l)th crossings is given by 



T, 



(fe+l/2) in 



exp 



(fe+l/2) 



<^mjm/4Vi (5i) 



or in the vector-matrix form 

c out _ j(k+l/2) c in 



(52) 



The change of the state across the fcth resonance is 
described by 



out _ Ak) in 



= A«< 



(53) 



(54) 



where for £,mj and £ ,m'j describing a pair of crossing 
levels ,, , is given by the LZ scattering matrix 

M, otherwise, if there is no crossing, it is A^ ,, , = 
S^^'Sm^mj. To formulate the condition of crossing, one 
has to consider the cases £' = ±1 separately. For £' = +1 
(descending line) the nuclear quantum number m'j = 

to ^ should be in the interval given by Eq. (|4"4"|) . for a 
crossing to be realized, then in the scattered state one 



has to/ 



,(-) 



m'j + k. For £' = — 1 (ascending line) 



the nuclear quantum number m'j — m) ^ should be in the 
interval given by Eq. (|45|) . for a crossing to be realized, 
then in the scattered state one has mj = rrij = m'j — k. 
Thus in the general case one obtains 



.4 



(55) 
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for £' = +1, max (—1 — k, —I) < m'j < min(7 — k, I), 



A (k) 
£,mi;i' ,m'i 



^.£'(^£,mf,m;+lt) 

)(56) 



for £ = — 1, max (—7 + fc, —7) < m' 7 < min(J+fc, I), and 



i(fe) 



'5 rn ' imi in all other cases. 



The evolution of the state across the whole grid is given 

by 

c out _ A (2/) T (2/-1/2) A (2/-1) , . . 

xA (-2/+l) T (-2J+l/2) A (-2/) c m j (57) 

where c ln is the state of the system before the leftmost 
crossing. The calculation can be made recurrent intro- 
ducing the scattering matrix after the fcth crossing: 



= L«< 



(58) 



where c ln is the initial state, same as in Eq. (|57|) . One 

has L(" 2/ ) = A^ 2/ ) 



-2/+1, 



L W =A (*) T (M/J) L (*-1) ) (59) 

-21 + 1, . . . , 21. In the final state it is 

L (2/ )c in , (60) 



-21 



1L.( 2Z ) being the full scattering matrix of the grid. 

If in the initial state the nuclear spin is in its ground 
state, application of Eq. (|57j) and calculation of the stay- 
ing probability yields the result of Eq. {J), no effect of 
nuclear spins. For 7=1/2 the problem can be solved 
analytically. If the nuclear spin is in its excited state, 
then the result is 



P = e 



+ exp 
1 — exp 



sA 2 
~F 2 

F 2 



1 — exp 



eA 2 
' F 2 " 



4 sin" 



eF 2 

7TA 2 



(61) 



In the fast-sweep limit e <C 1 this formula simplifies to 



P ^ e- £ + e z 



AK 
F 2 



4 sin 



eF 2 

7TA 2 



(62) 



i.e., at the linear order in e Eq. (fTJ) is robust. For 
A <C A x the envelope of Eq. (|6"T|) has a local maximum at 
e ~ A 2 1 A 2 with P max ~ 1. The asymptotic e ^> 1 behav- 
ior Eq. (|61[) is determined by the factor exp (— eA 2 / F 2 ^ 
that is slowly decaying for 4 < A x . The factor A 2 / F 2 
in the exponential is the square of the antidiagonal ele- 
ments of the matrix in Eq. (|37[) and it corresponds to the 
splitting that is much smaller than A. To the contrary, 
for A ^S> A x the result is close to P — er e . This means 
that the influence of nuclear spins strongly coupled to the 
electronic spin is small. 

In fact, analytical expressions, although too cumber- 
some, can be obtained for any J. The common feature 




25 £ 30 



FIG. 3: Asymptotic staying probability P for an electronic 
spin coupled to one nuclear spin / = 1/2, initially in its ex- 
cited state. The influence of strongly coupled nuclear spins 
becomes small, the curve A x /A = 1/3. 




FIG. 4: P for one nuclear spin I = 5/2, with A x /A = 1/3. 
Effect of the nuclear spin is large even for a small ratio A x /A 
due to large /. Averaging over the initial thermal distribution 
of nuclear spins (T — oo) reduces the effect. The incoherent 
approximation is defined in the next section. 



of the solutions for all I is a standard decay of P(e) for 
e ^> 1 in the case A ^> A x and an extremely slow decay in 
the case A <C A x . This becomes clear from the analysis of 
the LZ grids in both limiting cases, shown in Fig. [T] For 
A 3> A x the large splittings in Fig. [1] (a) line up horizon- 
tally, and in the slow-sweep limit the system cannot cross 
the dotted line moving along the ascending levels from 
left to right. It follows the exact levels and adiabatically 
turns down to the descending levels that leads to a full 
LZ transition, P « 0. On the contrary, for A <§; A x the 
large splittings in Fig. [1] (b) line up vertically. Although 
at slow sweep transitions at these crossings are adiabatic, 
there are smaller crossings at W < and W > where 
transitions are non-adiabatic, so that the electronic spin 
can end up in its both states. As a result, P(e) van- 
ishes only at extremely slow sweep, e ^> 1, at which the 
smallest crossings become adiabatic. 
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Next we present figures illustrating the dependence of 
the asymptotic staying probability P on the sweep-rate 
parameter e for the model with one nuclear spin. Fig. [3] 
shows the results for 7=1/2 and different ratios between 
A x and the electron-nuclear coupling A. The effect of 
nuclear spins is large for > A, in accordance with Eq. 

(ED- 

With increasing I the number of crossings in the LZ 
grid and thus the effect of the coupling to the nuclear spin 
on P increases, as one can see in Fig. 2J On the other 
hand, the effect is maximal for the nuclear spin initially 
in the mostly excited state and decreases in the case of 
the thermal distribution (T = oo) in the initial state. For 
I = 5/2 the effect is large even for A x /A = 1/3. This sug- 
gests that A x /A is not a proper parameter to describe rel- 
ative strength of different terms in the Hamiltonian. As 
commented upon below Eq. (|40p . the ratio between the 
sub-anti-diagonal to the dominant anti-diagonal terms is 



m,m+l 



A x /A ~ IA X /A 



(63) 



that is close to 1 in Fig. 2] Still in Mni2 the hyperfme cou- 
pling A is very strong, see comment below Eq. (|20[) . Ex- 
ternal transverse field of 2 T creates A x /A ~ 2/40 = 0.05, 
so that with I — 5/2 one has IA X /A ~ 0.12. Since this 
parameter enters squared the effective sweep parameter 
e mm i at sub-primary crossings, the direct effect of Mn 55 
spins on the LZ transitions in Mni2 should be small for 
not too strong magnetic fields. The physical reason for 
this is that nuclear spins do not have a sufficient dynamics 
to undergo transitions together with the electronic spin. 
The only thing that they can do is to create a static bias 
on the electronic spin, random because of the thermal 
distribution of nuclear spins. The latter becomes very 
important in combination with the dipole-dipole interac- 
tion between the electronic spins.— 



V. MANY NUCLEAR SPINS AND 
INCOHERENT APPROXIMATION 

Let us now consider the model of many nuclear spins 
coupled to the electronic spin. The simplest way to gen- 
eralize the results of the preceding section is to assume 
the equivalence of nuclear spins and the same coupling 
to the central electronic spin. In this case Eq. Ij20| is 
replaced by 



H = -\w{t)a z 



1 



-Act, 



(Aa z + A Z )I U 



A m L, 



Y x i tot 1 x 1 

(64) 

where Itot = Si=i * s the total nuclear spin. In fact, 
12 nuclear spins in Mni2 couple to the atomic spins with 
somewhat different coupling constants^ but this will be 
ignored for a moment. The length of the total nuclear 



H, (Itot)' 



0. Thus 



spin is dynamically conserved, 

the states of the nuclear subsystem can be classified by 
the value of the total spin / t ot < NI and its projection 
mi on some axis. This reduces the problem to the one 



0.008 



0.006 



0.004 



0.002 




• 7=5/2, JV= 12 
o /= 1/2, N= 146 
Gaussian 



0.000 



FIG. 5: Normalized distribution function for the total nuclear 
spin /tot- 



studied in the preceding section. In particular, for all 
nuclear spins in the ground state the nuclear subsystem 
will evolve as a single spin / to t = NI. In the realistic 
case of thermal distribution of nuclear spins one has to 
average the solutions for particular values of /tot and mj 
in the initial state over the distrubution of / to t and To/. 

The number of realizations 7V(/tot) of /tot can be 
computed recurrently^ If the total spin of a system 
of N nuclei is /tot, the total spin of its subsystem of 



N — 1 nuclei I' tot assumes the values |/t, 



min{/ to t + /, (N — 1)/}. Thus for the number of realiza- 
tions jV(Jtot ■ N) one can write 



min{7 tot +7,(7V-l)7} 



AA(/ t ot,A0 = 



E 



JVKt,JV-l). (65) 



The initial condition for this recurrence relation is 
Af(lU,2) = 1 for < / t ' ot < 21. The quantity Af(I to t) 
obeys the normalization condition 



NI 



J2 (2/tot + l)AA(/ tot ) = (2/+l) A '. 



(66) 



I tot =frac(iVI) 



For NI > 1, the quantity (2/ tot + l)AA(/ tot )/(2/+ 1) N is 
the high-temperature distribution function of the magni- 
tude of /tot and it is well approximated by 47r/ t 2 ot -F(Itot), 
where .F(Itot) is a normalized Gaussian function with re- 
spect to the three components of Itot<^ Thus for the 
asymptotic NI ^S> 1 form of the distribution function of 
/ to t normalized by 1 one has 



27T/ t , 



AA(/ to t) 

{2I+1) N (27TCT/) 3 / 2 



exp 



1 tot 
2(7 f 



(67) 



where 07 = (N/3)I(I + 1). It has a maximum at /tot = 
^/ctJ which is about 6 for Mni 2 (/ = 5/2, N = 12). 
Fig. [5]shows an agreement between the exactly computed 
M (Itat) / (21 + 1) N and its Gaussian approximation for 
1 = 5/2 and N = 12. This agreement improves for higher 
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7 = 5/2, N = \2 
or 

/= 1/2, N = 146 
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FIG. 6: P(e) for the electronic spin coupled to iV = 12 nuclear 
spins I = 5/2 (initially in a thermal state with T = oo) for 
different ratios A x /A and A z = 0. Results of the incoherent 
approximation are shown by solid lines of same color. 
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FIG. 7: P(e) for the electronic spin coupled to = 12 nuclear 
spins I = 5/2 for large ratios A x /A and A z = within the 
incoherent approximation. 



values of NI. One can see that practically the same result 
for Af(I tot )/(2I + 1) N can be achieved with I = 1/2 and 
N = 146. Thus in the model described by Eq. ([M]). the 
effect of 12 Mn 55 spins is the same as that of about 150 
protons. 

Summation over different values of itot makes calcula- 
tions lengthier. In addition, there are high values of itot 
with a small statistical weight that consume a lot of time 
but do not change the result. Thus one has to introduce 
a cut-off on itot- The results for I = 5/2 and N = 12 (or 
for 1 = 1/2 and N = 146) and A z = are shown in Fig. 
El Since the most probable value of itot m this case is 
between 5 and 6, deviations from the standard LZ effect 
are noticeable for A x /A down to 0.1 (see discussion at 
the end of preceding section) . 

One can see in Fig. [6] that summations over both itot 
and mj in the initial state lead to decreasing of quantum- 
mechanical oscillations in P(e). However, averaging out 
the oscillations is incomplete because for all values of itot 



the W intervals in the LZ grid are the same, although the 
number of crossings change. In real situations the oscilla- 
tions should be averaged out, mainly because of different 
couplings A to different nuclear spins that increases the 
number of crossings from (2i to t + I) 2 to (21 + 1) 2N . To 
take the phase averaging into account without increas- 
ing the size of the problem, one can use the incoherent 
approximation. In the case of well-separated resonances 
one can consider the occupation numbers of the states 



(68) 



The change of the state at the fcth resonance is described 
by Eq. ([55)1 . Then for the probabilities one has 



OUt _ A 



(fc) Ak)* in in* 

f ,mi ;£',m'j ,m, ,rn'/ L f ' ,m'j ,m'( ' 



(69) 



Neglecting interference effects or, in other words, averag- 
ing over phases amounts to the aproximation 



in „in* , „ c c 



(70) 



Then the change of the occupation numbers across a res- 
onance is described by 



i J i,m! — J -'£,mi;£',m' I * J £\m' I i 



where 



B 



(fc) 

£iirsi;f',m_. 



.4 



(fc) 

C,mi;£' ,m'j 



(71) 



(72) 



Incoherent approximation requires only a slight change 
of the computational method. Instead of Eq. ([5T)l one 
has 



(2/) 



(2/-1) 



»(-2/+l) 1 ,(-27) 1 



(73) 



In particular, for one nuclear spin 1=1/2 one obtains 
Eq. <jnU) with sin 2 (.. .) => 1/2. 

Results of the incoherent approximation are shown in 
Figs. [5] and [_o] by smooth curves. The results similar to 
those in Fig.[jo]hut for large ratios A x /A are shown in Fig. 
[7] within the incoherent approximation. For very large 
A x /A, the sub-primary splittings [left and right from the 
primary splittings at W — in Fig.[T](b)] are very small, 
so that for not too slow sweep transitions occur at W — 
and they become adiabatic with P w for e > 1. How- 
ever, with further increase of e transitions at crossings 
with sub-primary splittings begin, and the curve P(s) 
goes up again. Fig. [5] gives an idea of the dependence 
P(e) for a large nuclear spin. The apparent overall depen- 
dence is P(e) ~ C\ — Ci loge, seen in Fig. [8] over several 
decades in e. There are slow oscillations for A x /A = 100 
that tend to disappear for large /. 



VI. DISCUSSION 

The model of a two-level system coupled to one or 
many environmental spins considered above is one of ba- 
sic models illustrating many-body Landau-Zener effect 
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FIG. 8: P(e) for the electronic spin coupled to one large "nu- 
clear spin" for different ratios A x /A and A z = within the 
incoherent approximation. 

with many avoided level crossings forming a Landau- 
Zener grid. Exact results of the multilevel LZ effect such 
as the no-go theorem and DOBE formula apply to this 
model. In particular, if the environmental spins are ini- 
tially in their ground states, the asymptotic staying prob- 
ability is unchanged and given by Eq. ([1]). Accordingly 
no isotope effect exists in this case. However, such an 
effect occurs when preparing nuclear spins e.g., in a ther- 
mal ensemble. 

The model under consideration is relevant for spin tun- 
neling in molecular magnets, with coupling to nuclear 
spins, although comparison with experimental results re- 
quires taking into account the dipole-dipole interaction 
between electronic spins at the same time. The latter is 
more difficult and up to now could be done only in the 
fast-sweep limit £«liS 

Surprisingly, nuclear spins strongly coupled to the elec- 
tronic spin, such as 12 Mn 55 nuclear spins I = 5/2 in 
Mni2, do not significantly change the asymptotic sur- 
vival probability P(e) for the electronic spin. This fol- 
lows from the analysis of the Landau-Zener grid, Fig. [1] 
(a). In physical terms, nuclear spins in the case A x <^ A 
do not have a sufficient dynamics to undergo transitions 
together with the electronic spin. They rather act on 
the latter as static random fields that have no effect by 
themselves but do significantly reduce the influence of 
the DDI on the LZ effect^ 

To the contrast, protons in molecular magnets can in- 
fluence the LZ effect directly. Although the nuclear mag- 
neton fi n is small, the energy of the dipole-dipole inter- 



action with nuclear spins can exceed tunnel splitting A. 
For instance, the ground-state splitting in Fes i n zero 
field is only A/ fee — 10~ 7 K. The energy bias on the 
electron spin from a proton at distance r is given by 
W ~ 2Sgn B fi p /r 3 , where g = 2 and fi p = 2.79/i„. One 
has A ~ W at the distance r& — 72 A. The volume of 
a corresponding sphere comprises about 400 unit cells of 
a molecular magnet. As each molecule, occupying a unit 
cell, comprises about 120 protons^ the number of pro- 
tons creating bias W > A is about N = 50000. These 
protons are typically weakly coupled to the electronic 
spin (^4 <C Aa;) since their Zeeman interaction with the 
external field greatly exceeds their NDDI with electronic 
spins. Applying an external magnetic field, one can dras- 
tically increase A and thus reduce N. 

Effect of protons in principle could be tackled with the 
method described in Sec. [II] However, as each proton 
couples to the electronic spin with its own coupling con- 
stant A, the number of lines in the LZ grid, 2 , is too 
large to implement a working computational algorithm. 
The first expected effect of so many crossings is complete 
averaging out of quantum-mechanical phase oscillations 
that are still seen in Fig. [51 This effect can easily be 
accounted for by the incoherent approximation of Sec. 
fVl After that different couplings A become less impor- 
tant, and one can get a qualitative idea of the effect of 
protons from the analysis in the preceding section. The 
most probable combined nuclear spin / to t in the distri- 
bution of Eq. (|67p is y/aj w 110 that is of the same 
order as the large "nuclear spins" simulated in Fig. [5] 
The resulting P{e) is decaying extremely slowly because 
of a strong non-adiabaticity at large e induced by weak 
avoided crossings outside the ellipse in Fig. As the 
width of the distribution of /tot is of the same order as 
the most probable value of I tot , slow oscillations in Fig. 
[2] should be averaged out after summation over all /tot ■ 
One can see that the results for P(e) in the case A <C A. x 
slowly approach 1/2. This is in accord with the results of 
Ref. |46| where the influence of environment was modeled 
by a density-matrix equation with dephasing. 

It does not make sense, however, to further elaborate 
on the effect of protons in this paper because in the 
practical case 4 < A x coupling to electronic spins only 
slightly perturbs protons. This might be an indication of 
a possibility to solve the problem by another and more 
efficient method. 

The work by D. A. G. was supported by the NSF grant 
No. DMR-0703639. 
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